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Abstract 

Using the entire CLEO-c tp (3770) — > DD event sample, corresponding to an integrated lumi- 
nosity of 818 pb~^ and approximately 5.4 million DD events, we present a study of the decays 
D^ — 7r~e+fe, D^ — K'e^Uf,, D^ — 7r°e+fe, and D^ — K^e^Ve- Via a tagged analysis tech- 
nique, in which one D is fully reconstructed in a hadronic mode, partial rates for semileptonic 
decays by the other D are measured in several bins. We fit these rates using several form factor 
parametcrizations and report the results, including form factor shape parameters and the branching 
fractions ^ (£»o ^ TT'e+Ve) = (0.288±0.008±0.003)%, B {D^ K-e+Ue) = (3.50±0.03±0.04)%, 
B {D+ 7r°e+z/e) = (0.405 ± 0.016 ± 0.009)%, and B {D+ K^e+v^) = (8.83 ± 0.10 ± 0.20)%, 
where the first uncertainties are statistical and the second are systematic. Taking input from lat- 
tice quantum chromodynamics (LQCD), we also find \Vcd\ = 0.234 ± 0.007 ± 0.002 ± 0.025 and 
\Vcs\ = 0.985 ± 0.009 ± 0.006 ± 0.103, where the third uncertainties are from LQCD. 



* Deceased 
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I. INTRODUCTION 



Semileptonic decays are an excellent environment for precision measurements of the 
Cabibbo-Kobayashi-Maskawa (CKM) [l|, 0] matrix elements. However, because these decays 
are governed by both the weak and strong forces, extraction of the weak CKM parameters 
requires knowledge of strong interaction effects. These can be parameterized by form fac- 
tors. Techniques such as lattice quantum chromodynamics (LQCD) offer increasingly precise 
calculations of these form factors, but as the uncertainties in the predictions shrink, experi- 
mental validation of the results becomes increasingly important. Because the magnitudes of 
CKM matrix elements \ Vcd\ and \ Vcs\ are tightly constrained by CKM unitarity, semileptonic 
decays of D mesons provide an excellent testing ground for the new theoretical predictions. 
The relevance of tests using charm decays is increased by the similarity of D meson decays 
to those of B mesons, where QCD calculations are critical to extractions of \Vub\ in exclusive 
B semileptonic decays j3|. 

We present a study of the decays D° — > n~e~^h'e, K~e^Uf,, — > n^e^Ve, and — > 

K^e^Ue (with charged conjugate modes implied throughout this article) using 818 pb~^ of 
?/'(3770) DD data collected by the CLEO-c detector. Taking advantage of the fact 
that D mesons produced near the '?/'(3770) resonance are produced solely as part of DD 
pairs, we follow a tagged technique pioneered by the Mark III Collaboration [3| and used in 
semileptonic analyses of smaller portions of CLEO-c data jl, 0, 0] . Hadronically decaying D 
tags are first reconstructed; one then looks for the D decays of interest in the remainder of 
each event. This strategy suppresses backgrounds and provides an absolute normalization 
for branching fraction measurements. 

For semileptonic decays such as those of interest here, in which the initial and final state 
hadrons are pseudoscalars and the lepton mass is negligibly small, the strong interaction 
dynamics can be described by a single form factor /+ (g^), where is the invariant mass of 
the lepton-neutrino system. The rate for a D semileptonic decay with final state meson P 
is given by 



Vcd{s) 



.p3 



where Gp is the Fermi constant, Vcd(s) is the relevant CKM matrix element, p is the mo- 
mentum of the daughter meson in the rest frame of the parent D, and X is a multiplicative 
factor due to isospin, equal to 1 for all modes except D~^ — >• yr^e+z/e, where it is 1/2. The 
primary measurements described here are the partial decay rates AF = / ^dq^ in seven q^ 
bins each for D^ -K^e^Vf. and D^ -K^e^Vf, and nine q^ bins each in D^ K^e^Ve and 
D^ K^e^Vf.- We fit the AF using several parameterizations of /+ (g^), extracting form 
factor shape parameters, measurements of Vcd{s) /+(0), and branching fractions. Taking 
estimates of /+(0) from theory, we also extract \Vcd\ and \Vcs\- 

The remainder of this article is organized as follows: the CLEO-c detector and event 
reconstruction are described in Sec. [Tll Measurements of partial rates and their system- 
atic uncertainties are detailed in Sees. IIIII and IIVI respectively. Extractions of branching 
fractions, form factor shapes, and CKM parameters are reviewed in Sec. El and Sec. I VI I 
summarizes our results. 
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II. DETECTOR AND EVENT RECONSTRUCTION 



The CLEO-c detector has been described in detail elsewhere @, i, [13 . The 53-layer track- 
ing system, composed of two drift chambers covering 93% of the solid angle and enclosed by 
a superconducting solenoid operating with a 1-Tesla magnetic field, provides measurements 
of charged particle momentum with a resolution of ~ 0.5% at 700 MeV/c. The tracking 
chambers also supply specific-ionization (dE/dx) information, which is combined with in- 
put from the Ring-Imaging Cerenkov (RICH) detector to provide excellent discrimination 
between charged pions and kaons. A 7784 crystal cesium-iodide calorimeter covering 95% 
of the solid angle provides photon energy resolution of 2.2% at E = 1 GeV, with a vr" mass 
resolution of about 6 MeV/c^, and contributes to positron identification. 

The entire CLEO-c 'ijj{3770) data sample has an integrated luminosity of 818 pb"^, equiv- 
alent to approximately 5.4 million DD events. The data were collected at center-of-mass 
energies near 3.774 GeV with a RMS spread in beam energy of approximately 2.1 MeV. 
Events collected at this energy, approximately 40 MeV above the DD production threshold, 
are composed primarily of D^D^, D~^D~ , and non-charm continuum final states. 

GEANT-based [ll| Monte Carlo (MC) simulations are used to determine reconstruction 
efficiencies, develop line shapes for yield extraction fits, and conduct tests of the analysis 



procedure. Final state radiation (FSR) is simulated using PHOTOS [12| version 2.15 
with the option to simulate FSR interference enabled. A sample of generic '?/'(3770) — *• 
DD events, generated using EvtGen ^] and corresponding to approximately 20 times the 
data luminosity, was generated using input from Ref. [13], combined with CLEO-c results 
using the initial 281 pb^^ data sample where appropriate. This sample, along with samples 
of simulated e+e" ^ = u,d, or s), e^e" t^t~, and e~^e~ ip{2S)j events 

corresponding to five times the data luminosity, is referred to as "generic MC" for the 
remainder of the article. We also use a sample of ?/'(3770) —>■ DD events in which the D 
meson decays to one of the four studied semileptonic modes and the D decays to one of the 
hadronic final states used in tag reconstruction. This sample is referred to as "signal MC." 
In both the generic and signal MC samples, the semileptonic decays are generated using the 
modified pole parameterization [is!] (see Sec. IV AD with parameters fixed to those measured 
in the initial 281 pb~^ of CLEO-c data [B], 1^. All simulations are corrected for small biases 
in the positron, charged hadron, and 7r° identification efficiencies. 

Charged pions and kaons are identified from drift chamber tracks with momentum greater 
than 50 MeV/c and with | cos^| < 0.93, where 6 is the angle between the track and the beam 
axis. Charged track reconstruction efficiencies are approximately 84% for kaons and 89% 
for pions; lost tracks within |cos^| < 0.93 are almost exclusively due to particle decay in 
fiight and material interaction in the drift chambers. Pions and kaons are distinguished 
using a combination of specific ionization measurements and, if the track momentum is 
greater than 700 MeV/c, RICH detector information. For all other tracks, hadron identity is 
determined using specific ionization information only. Given a properly reconstructed track, 
hadron identification efficiencies are approximately 95%, with misidentification rates of a few 
percent. Identical hadron selection criteria are used in tag and semileptonic reconstructions. 

Neutral pion candidates are reconstructed via vr^ —>■ 77. Photon candidates are identified 
from energy depositions in the calorimeter greater than 30 MeV using shower shape informa- 
tion. The invariant mass of each pair of photon candidates is calculated using a kinematic 
fit that assumes the photons originate at the center of the detector. This mass is required 
to be within three standard deviations (3cr) of the nominal 7r° mass, where a is determined 
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TABLE I: AE = E^t^g — -E'beam requirements for tag reconstruction. 



Mode 






Requirement (GeV) 


D° 


K+-K- 




|A£;| < 0.030 


D° 


K+ir- 




-0.050 <AE< 0.044 






TT TT^ 


\AE\ < 0.020 


D- 






|A£;| < 0.0232 


D- 






< 0.0276 


D- 


KItt- 




|A^| < 0.0272 


D- ^ 




^0 


A^ < 0.0366 


D- 






A^ < 0.0159 


D- 


K+K 




A^ < 0.0138 



from the kinematic fit. The resulting vr" energy and momentum from the fit are used in 
further event analysis. Efficiencies for tt^ reconstruction vary from 40% at a momentum of 
100 MeV/c to 60% at 900 MeV/c. If multiple neutral pion candidates are reconstructed 
opposite the tag, the candidate with the mass closest to the nominal vr*^ mass is chosen. 

Neutral kaon candidates are reconstructed via Kg —>■ tt~^7[~ using vertex-constrained fits 
to pairs of oppositely charged intersecting tracks. The invariant mass of the tt+tt" candidate 
is required to be within 12 MeV/c^ of the nominal Kg mass. This procedure results in a 
Kg mass resolution of 2 - 2.5 MeV/c^ and a Kg reconstruction efficiency of about 94%. If 
multiple Kg candidates are reconstructed opposite a tag, the candidate with mass closest 
to the nominal Kg mass is chosen. 

Positron candidates are identified from tracks with momentum greater than 200 MeV/c 
and within the solid angle |cos^| < 0.9. Positrons are selected using a combination of 
specific ionization, calorimetry, and RICH detector information. The efficiency for positron 
identification is about 50% at the low momentum threshold of 200 MeV/c, rises sharply 
to 92% at 300 MeV/c, and varies by a few percent as a function of momentum beyond 
300 MeV/c. Roughly 0.1% of charged hadrons satisfy the positron identification criteria. 
Positron momentum resolution is degraded by FSR. We reduce this effect by identifying 
bremsstrahlung photon candidates in the calorimeter within 5° of the positron candidate 
track and adding their 4-momenta to that of the positron candidate. Such photons must 
have energy greater than 30 MeV and no associated track reconstructed in the drift chamber. 

Tag candidates are reconstructed in three D° decay modes (Z)° K~^'k~, — > K^tt~7i^, 
and —>■ K^-K^-K^-K^) and six D~ decay modes {D^ K^tt^ti^, K^tt^tt^tt^, 
D- K^s^~, D- KItx-h^, D~ KI-K-TI-TI+, and D' K+K'ti-). Backgrounds 
are suppressed by requiring that AE = i^tag ~ -^bcam satisfy the requirements given in 
Table [H These cuts correspond to approximately ±4cr, with a depending on the de- 
cay mode. Backgrounds are further reduced using the beam-constrained mass, Mbc = 

1 /2 

(-^beam/c^ — |Ptag|^ /c^) , whcrc -Ebeam IS the beam energy and Ptag is the total measured 
momentum of the tag candidate. The Z)° tag candidates must satisfy 1.858 < Mbc < 1-874 
GeV/c^, while D~ tag candidates are required to have 1.8628 < Mbc < 1.8788 GeV/c^. 
In events with multiple tag candidates, the one candidate per mode and per D flavor with 
reconstructed energy closest to the beam energy is chosen. 

Semileptonic candidates are formed from positron and hadron candidate pairs. Although 
the semileptonic neutrino daughter is not detected, its energy and momentum can be inferred 
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from the missing energy £^miss and momentum Pmiss of the event: 

-E^miss -E^beam ^he (2) 

and 

Pmiss = ~Ptag ~ P/ie; (3) 

where the energy Ehe and momentum Phe of the hadron-positron system are constructed 
from the measured energy and momenta of the hadron, positron and any bremsstrahlung 
photon candidates. The tag momentum PJg^g is formed from the measured tag momentum 
with the magnitude constrained using the beam energy and D mass: P[g,g — [(-E'beam/c)^ — 

{cm o^Y^'^P tug- All momentum vectors are boosted to the center-of-mass frame by correcting 
for the small e+e" net momentum due to the beam crossing angle (~ 3 mrad). 

Semileptonic backgrounds are reduced by requiring that the variable U, defined as 

U = -Emiss c| Pmiss I, (4) 

satisfy —0.10 < U < 0.24 GeV for each candidate. Additionally, the positron and hadron 
from the semileptonic decay are required to have opposite charge in D° — > Tr~e'^i'e and 
K'e^Ve candidates, while the positron and tag are required to have opposite charge 
in D+ 7r°e+z/e and D+ -> K^e+u^. 

Semileptonic candidates are partitioned into several bins, where the reconstructed q'^ 
is determined based on measurements of the positron and neutrino: 



q^^l-(E^ + E,f-\\P^ + P,\\ (5) 



The neutrino energy is taken to be the missing energy of the event, while the neutrino 
momentum is equated with the missing momentum scaled so that |Pmiss| = -E'miss- Because 
-Smiss does not require measurements of the tag decay, it is a better measured quantity than 
I Pmiss I; constraining the neutrino momentum in this manner thus improves the resolution 
in q^. TT~e^Ve and T^^e^Ve candidates are divided into seven q^ bins, with 

boundaries defined by [0,0.3), [0.3,0.6), [0.6,0^9), [0.9,1.2), [1.2,1.5), [1.5,2.0), and [2.0, oo) 
GeV^/c^. In the D° — > K't^v^ and — > K^e^v^ modes, nine q^ bins are used, defined 
by [0,0.2), [0.2,0.4), [0.4,0.6), [0.6,0.8), [0.8,1.0), [1.0,1.2), [1.2,1.4), [1.4,1.6), and [1.6, oo) 
GeVVcl 



III. EXTRACTION OF PARTIAL RATES 

In order to measure the partial rates we first determine the number of observed tags 
A^tag^'" ("tag yields") in each tag mode a. This is related to the number of tags produced in 
mode a, A^^ag^ by 

A 7-obs,a 

^ug = (6) 

^tag 

where e^g^g is the reconstruction efficiency for tag mode a. We then determine the number 
of events with both a tag and semileptonic candidate. These "signal yields" are 
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determined separately for each tag mode a and bin j. The signal yields are related to 
the number of tag-semileptonic combinations produced in each bin, n°, by 

<"" = E<ei.' (7) 

i 

where e^- are the elements of a matrix that describes the efficiency and smearing across 
bins associated with tag and semileptonic reconstruction. As the number of tag-semileptonic 
combinations produced is a function of the number of tag decays and the differential semilep- 
tonic decay rate, dV /dq^, we can rewrite Eq. ([7]) as 

<^•" = iVtVz.E4/|^^^^ (8) 

where td is the D lifetime and the integration is over the width of bin i. Combining 
the above equations and solving for the differential rate, we obtain a simple formula for 
extracting the partial rates: 

Mr , 2 

J I (iq 

The small correlation between signal and tag yields that arises from the signal yields being 
a subset of the tag yields is neglected. 

The following sections describe the extraction of the tag yields iV°ag tagging efficiencies 
e"jjg, signal yields nf^^'°', and signal smearing and efficiency matrices e^-. With all of these 
numbers in hand, we then extract the partial rates. 



A. Tag Yields and EfBciencies 

The tag yields A^t°ag*'" are obtained separately for each mode a by performing fits to 
beam-constrained mass distributions. The fitting procedure has been described in detail in 



17l | and involves an unbinned likelihood maximization. True tag decays are modeled using a 
function specially designed to take into account the natural '?/'(3770) line shape, beam energy 
resolution, momentum resolution, and initial state radiation (ISR) effects. Tag backgrounds 



are modeled using an ARGUS function [18|, modified so that the power parameter is allowed 



to float [17||. The fits, shown in Fig. [T], are performed over a wide 1.83 < Mbc < 1-89 
GeV/c^ window. Tag yields, given in Table UTt are obtained by subtracting the backgrounds 
estimated by the fits from event counts in data inside the narrower Mbc signal regions. Also 
shown in Table [III are tagging efficiencies, which are obtained by fitting generic MC Mbc 
distributions with the same procedure used in data. 



B. Signal Yields and Efficiencies 

Signal yields are extracted from distributions of U, defined in Eq. (jl]). Events in which 
both the tag and semileptonic decay have been correctly reconstructed, leaving only an 
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3070109-001 




1.84 1.86 1.88 1.84 1.86 1.88 1.84 1.86 



(GeV/c2) 

FIG. 1: Mbc distributions in data (points), with fits (solid lines) and background contributions to 
fits (dotted lines). The vertical lines show the limits of the Mbc signal regions. 

undetected neutrino, are expected to peak at U = 0, with the shape of the distribution being 
approximately Gaussian due to detector resolution. Misreconstructed events and background 
modes generally have non-zero U values. Properly reconstructed decays are separated from 
backgrounds using an unbinned maximum-likelihood fit, executed independently for each 
semileptonic mode, each tag mode, and each bin. A sample of the U distributions is 
shown in Figs. [2] - [51 

For the fit, the U distribution of signal candidates is taken from signal MC samples. 
While the U resolution in data is approximately 12 MeV for the modes with only charged 
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TABLE II: Tag yields and statistical uncertainties in data and tag reconstruction efficiencies. 



Mode 






Yield 


Efficiency(%) 


1)0 ^ 


K+TT- 




149616 ± 392 


65.32 


D° 


K+TT- 




284617 ± 589 


35.15 


D^^ 


K+TT- 




227536 ± 517 


45.55 


D- - 


K+TT- 




233670 ± 497 


55.42 


D- - 


K+TT~ 




69798 ± 330 


27.39 


D- - 


K^TT- 




33870 ± 194 


51.10 


D- - 




vrO 


74842 ± 357 


28.74 


D 






49117 ±323 


43.58 


D- 


K+K 




19926 ± 171 


42.07 



tracks in the final state {D*^ vr^e+i/g, K~e'^i'e, and K^e^Ve) and 25 MeV for 

— >■ TT^e^Vei tlie distributions are sliglitly narrower in MC simulation. To accommodate 
the poorer U resolutions in data, the MC distributions are convolved with a double Gaus- 
sian with parameters fixed for each semileptonic mode to the values that maximize the fit 
likelihoods summed over all bins and tag modes. The smearing functions are dominated 
by central Gaussians with widths of approximately 6 MeV in — > K~e^Vf., —>■ 7i~e~^i'e, 
and — > K^e'^Ug and 13 MeV in — > 7r^e~^i>e. The secondary Gaussians have normal- 
izations of 3% - 7% of the central Gaussian and have widths of 30 - 35 MeV. The overall 
normalization of the corrected signal distribution is allowed to float in each fit. 

The background distributions used in the fits are taken from generic MC samples. Back- 
grounds arise from both non-DD and DD events. The —>■ TT~e^Ve and — >■ ir^e^Uf. 
modes are subject to large backgrounds from K~e^Ve and K^e^Vei respec- 

tively. To allow for variations in vr" and vr" fake rates between the data and MC simulations, 
the normalizations of these components are fixed to the values that minimize the fit likeli- 
hood summed over all bins and tag modes. The — >■ ir'e^Ve mode is also subject to a 
large background from p~e~^i>e events; the normalization of this background is fixed 

using the known branching fraction and the tag yields in data and MC samples. The remain- 
ing DD backgrounds occur due to misreconstruction of either the semileptonic or tag decay, 
although the largest backgrounds are composed of events with a correctly reconstructed tag 
but misreconstructed semileptonic decay. In each semileptonic mode, all of the DD back- 
grounds not discussed above are combined into a single background distribution with the 
normalization allowed to float in each fit. The normalization of the non-DD distribution is 
fixed using the ratio of luminosities in data and MC samples. 

Because each bin and tag mode is treated separately, the total numbers of fits for 
vr-e+z/e, D° R-e+Ue, D+ n°e+Ue, and D+ K'^e+Ue are 21, 27, 42, and 54, 
respectively. Figs. [2] - [5] show four individual fits, as well as the summed fit results for all 

bins and tag modes. 

The signal efficiency matrix elements eg-, as defined in Eq. ([7]), are obtained from signal 
MC simulations and corrected for previously determined deviations from positron, charged 
hadron, and vr" identification efficiencies in data. Each eg- gives the fraction of events gener- 
ated in bin j with tag mode a that are reconstructed in g^ bin i with the same tag. The 
efficiency matrix thus accounts for reconstruction of both the signal and tag decays. The 



D~^ Ti^e^Ve efficiencies include the vr*^ 77 branching fraction 19(] and D^ — >- K^e 
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FIG. 2: U distributions in data (points) for — > Tr~e'^U(., with fit results (liistograms) showing 
signal (clear) and background components: p~e^Ue (darkest gray), — > K~e^Ve (lightest 

gray), other DD (medium gray), and non-DD (black). 
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FIG. 3: U distributions in data (points) for — > K e^Ve, with fit results (histograms) showing 
signal (clear) and background components: DD (gray) and non-DD (black). 



efficiencies include the K'g fraction of the K*^ and tt+tt^ branching fraction [19]. In 

total, there are eighteen efficiency matrices - one for each tag and semileptonic mode combi- 
nation. Tables IIIII and [IV] provide four examples of these matrices. The diagonal elements, 
giving the efficiency for the tag and semileptonic decays to be reconstructed in the correct 
bin, vary from 5% - 50% depending on semileptonic mode, tag mode, and q^. The neigh- 
boring off-diagonal elements, giving the efficiencies for the tag and semileptonic decay to 
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FIG. 4: U distributions in data (points) for — > with fit results (histograms) showing 

signal (clear) and background components: — > K^e'^Ue (light gray), other DD (dark gray), 
and non-DD (black). 
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FIG. 5: U distributions in data (points) for K^e'^Ve, with fit results (histograms) showing 

signal (clear) and background components: DD (gray) and non-DZ) (black). 



be reconstructed in the wrong bin, range between 1% and 10% of the diagonal elements. 
The signal yields summed over tag modes both before and after correction by these matrices 
are shown in Table El 
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TABLE III: Selected efficiency matrices e^j in percent for n^e^Uf. and K~e^Ve- 

The column gives the true bin while the row gives the reconstructed ("Rec") bin i. The 

elements account for the reconstruction efficiencies of both the tag and the semileptonic decay. 
The statistical uncertainties in the least significant digits are given in the parentheses. 

TT- e+i^e , — >■ K+-K- 

Rec g2 True (GeVVc") 

(GeV^/c"^) [0,0.3) [0.3,0.6) [0.6,0.9) [0.9,1.2) [1.2,1.5) [1.5,2.0) [2.0,oo) 

[0, 0.3) 41.25(34) 1.19(8) 0.02(1) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 

[0.3, 0.6) 0.76(6) 42.57(36) 1.55(10) 0.01(1) 0.00(0) 0.00(0) 0.00(0) 

[0.6, 0.9) 0.04(1) 1.12(8) 44.65(38) 1.54(10) 0.02(1) 0.00(0) 0.00(0) 

[0.9, 1.2) 0.02(1) 0.08(2) 1.09(8) 44.77(41) 1.37(10) 0.03(1) 0.00(0) 

[1.2, 1.5) 0.01(1) 0.03(1) 0.09(2) 1.33(9) 46.11(44) 0.91(8) 0.00(0) 

[1.5, 2.0) 0.01(1) 0.02(1) 0.02(1) 0.11(3) 1.20(10) 47.01(40) 0.74(8) 

[2.0, oo) 0.00(0) 0.00(0) 0.01(1) 0.02(1) 0.04(2) 0.56(6) 47.28(48) 

r>° — *■ K-e+u^,D° K+-K--K° 

Rec g2 True (GeVVc*) 

(GeV^/c"t) [0,0.2) [0.2,0.4) [0.4,0.6) [0.6,0.8) [0.8,1.0) [1.0,1.2) [1.2,1.4) [1.4,1.6) [1.6, oo) 

[0, 0.2) 19.70(4) 0.80(1) 0.03(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 

[0.2, 0.4) 0.45(1) 19.80(5) 1.03(1) 0.03(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 

[0.4, 0.6) 0.01(0) 0.54(1) 20.58(5) 1.12(1) 0.03(0) 0.01(0) 0.00(0) 0.00(0) 0.00(0) 

[0.6, 0.8) 0.01(0) 0.02(0) 0.61(1) 21.32(6) 1.12(2) 0.03(0) 0.01(0) 0.00(0) 0.01(0) 

[0.8, 1.0) 0.01(0) 0.01(0) 0.03(0) 0.63(1) 21.92(6) 1.03(2) 0.01(0) 0.00(0) 0.00(0) 

[1.0, 1.2) 0.00(0) 0.01(0) 0.02(0) 0.03(0) 0.59(1) 21.64(7) 0.95(2) 0.01(0) 0.00(0) 

[1.2, 1.4) 0.00(0) 0.00(0) 0.01(0) 0.01(0) 0.01(0) 0.51(1) 21.08(9) 0.88(3) 0.01(0) 

[1.4, 1.6) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.01(0) 0.39(1) 20.05(11) 0.79(4) 

[1.6, oo) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.25(1) 16.72(17) 



C. Partial Rate Results 

Using the tag yields, tag efficiencies, signal yields, and signal efficiency matrices, we solve 
Eg . (p) for the partial rates in each bin and tag mode, AF". The procedure detailed 
in |20| is used to calculate uncertainties and correlations in the inverted efficiency matrices. 
We then average the resulting AT" over tag modes, obtaining AFj. Statistical covariance 
matrices detailing the uncertainties on the AFj are also calculated and are available in the 
Appendix. Within each semileptonic mode, there are small correlations across bins that 
arise from the smearing in q^. Both the individual and tag-averaged partial rates are shown 
in Fig. O The tag-averaged partial rates are also given in Table El 



D. Tests of Partial Rate Results 

Calculating the partial rates separately for each tag mode allows for a test of the con- 
sistency of results across different tag modes. To quantify the tag mode agreement, we 
calculate a for each semileptonic mode: 

xW^EE^^^f^. (10) 

where erf is the statistical uncertainty on AF". This quantity is expected to have a 
distribution, with mean ridof and variance 2?2dof, where the number of degrees of freedom 
is given by ndof = ''^g2 ^ms x ('^tag modes ~ !)• Table IVT] gives the measured the number 
of degrees of freedom, and the probability for each mode. These values show that the 
semileptonic rates agree well across tag modes. 
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TABLE IV: Selected efficiency matrices eij in percent for — > n^e'^Uf, and — > K^e ' v^. 
Tlie column gives the true bin while the row gives the reconstructed ("Rec") bin i. The 
elements account for the reconstruction efficiencies of both the tag and the semileptonic decay. 
The statistical uncertainties in the least significant digits are given in the parentheses. 



Rec (jf^ 
(GcVVc'') 



True g2 (GeVVc*) 
[0,0.3) [0.3,0.6) 



[0.6,0.9) [0.9,1.2) [1.2,1.5) [1.5,2.0) [2.0, oo) 



[0, 0.3) 



[0.3, 
[0.6, 
[0.9, 
[1.2, 
[1.5, 
[2.0, 



0.6) 
0.9) 
1.2) 
1.5) 
2.0) 
oo) 



22.44(20) 
1.23(5) 
0.03(1) 
0.02(1) 
0.02(1) 
0.02(1) 
0.17(2) 



0.83(5) 
21.69(21) 
1.62(6) 
0.03(1) 
0.03(1) 
0.03(1) 
0.19(2) 



02(1) 

02(5) 

23(22) 

75(7) 

06(1) 

04(1) 

31(3) 



0.00(0) 
0.01(1) 
1.14(6) 
21.12(23) 
1.61(7) 
0.13(2) 
0.47(4) 



0.00(0) 
0.00(0) 
0.01(1) 
1.05(6) 
19.72(25) 
1.47(7) 
0.70(5) 



0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.65(4) 
20.50(22) 
1.65(7) 



0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.49(5) 
22.81(27) 



Rec 
(GeVVc^) 



[0, 0.2) 



[0.2, 
[0.4, 
[0.6, 
[0.8, 
[1.0, 
[1.2, 
[1.4, 
[1.6, 



0.4) 
0.6) 
0.8) 
1.0) 
1.2) 
1.4) 
1.6) 
oo) 



D+ 

True (j2 (GeVVc*) 

[0,0.2) [0.2,0.4) [0.4,0.6) [0.6,0.8) [0.8,1.0) [1.0,1.2) [1.2,1.4) [1.4,1.6) [1.6, oo) 
5.06(3) 
0.11(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 



0.21(1) 
4.98(3) 
0.15(1) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 



00(0) 
24(1) 
09(3) 
16(1) 
00(0) 
00(0) 
00(0) 
00(0) 
01(0) 



0.00(0) 
0.00(0) 
0.25(1) 
5.12(3) 
0.15(1) 
0.01(0) 
0.00(0) 
0.00(0) 
0.01(0) 



0.00(0) 
0.00(0) 
0.01(0) 
0.28(1) 
5.13(4) 
0.15(1) 
0.01(0) 
0.00(0) 
0.01(0) 



0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.26(1) 
5.14(4) 
0.12(1) 
0.01(0) 
0.01(0) 



0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.24(1) 
5.29(5) 
0.12(1) 
0.01(0) 



0.00(0 
0.00(0 
0.00(0 
0.00(0 
0.00(0 
0.00(0 
0.22(1 
5.34(7 
0.08(1 



0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.00(0) 
0.01(0) 
0.26(3) 
5.45(11) 



TABLE V: Signal yields, both raw (n° '^'") and corrected for q'^ smearing and reconstruction 
efficiency (nf), and partial rates (AFj). Statistical uncertainties in the least significant digits are 
given in parentheses. 

g2 (GeVVc") [0,0.3) [0.3,0.6) [0.6,0.9) [0.9,1.2) [1.2,1.5) [1.5,2.0) [2.0, oo) 

Raw yield 251(17) 232(16) 204(15) 194(15) 161(13) 173(14) 159(13) 

Corrected yield 858(60) 795(59) 636(51) 612(50) 495(45) 532(46) 505(45) 

Partial rate (ns~^) 1.39(10) 1.22(9) 1.02(8) 0.98(8) 0.79(7) 0.84(7) 0.80(7) 

g2(GeVVc'') [0,0.2) [0.2,0.4) [0.4,0.6) [0.6,0.8) [0.8,1.0) [1.0,1.2) [1.2,1.4) [1.4,1.6) [1.6, oo) 

Raw yield 2751(54) 2541(51) 2325(49) 2034(46) 1662(42) 1243(36) 904(31) 496(23) 167(13) 

Corrected yield 11290(233) 10026(222) 8629(203) 7432(186) 5867(163) 4485(143) 3393(125) 1983(98) 781(68) 

Partial rate (ns-i) 17.82(36) 15.83(35) 13.91(32) 11.69(29) 9.36(26) 7.08(22) 5.34(19) 3.09(15) 1.28(11) 

D+ 7r°e+i/e 

9^ (GeV^/c*) [0,0.3) [0.3,0.6) [0.6,0.9) [0.9,1.2) [1.2,1.5) [1.5,2.0) [2.0, oo) 

Raw yield 148(13) 141(13) 124(12) 122(12) 100(11) 107(12) 96(13) 

Corrected yield 799(80) 748(82) 665(80) 640(77) 570(80) 625(82) 460(86) 

Partial rate (ns-^) 0.71(7) 0.66(7) 0.56(6) 0.57(6) 0.48(6) 0.54(7) 0.37(7) 

D+ K.°e+Ue 

q^{GeV^/c-^) [0,0.2) [0.2,0.4) [0.4,0.6) [0.6,0.8) [0.8,1.0) [1.0,1.2) [1.2,1.4) [1.4,1.6) [1.6, oo) 

Raw yield 1704(42) 1511(40) 1389(38) 1229(36) 912(31) 809(29) 514(24) 275(17) 124(12) 

Corrected yield 10090(282) 8732(271) 7934(256) 6951(240) 5101(207) 4511(190) 2812(152) 1412(106) 625(73) 

Partial rate (ns-i) 17.79(47) 15.62(45) 14.02(43) 12.28(40) 8.92(34) 8.17(32) 4.96(25) 2.67(18) 1.19(13) 



As a test of the signal yield fits, we compare the observed and predicted distributions in 
three variables: the lepton momentum, the hadron momentum, and the angle of the virtual 
W in the D rest frame relative to the positron in the W rest frame. Fig. [7] shows these 
distributions for data and MC candidates with \U\ < 60 MeV, with signal and background 
MC distributions normalized using the U fits. All of the distributions show good agreement 
between data and MC simulations. 
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FIG. 6: Partial rates for each semileptonic mode. The points show measurements in each tag 
mode; the histograms show the partial rates averaged over all tag modes. 



TABLE VI: of partial rates across tag modes, with number of degrees of freedom and x 
probability. 



Semileptonic 


mode 




^dof Pix^) 


D° ^ TT-e 




12 


14 61% 


L>° ^ K-e 




21 


18 28% 


D+ -K^e 


'^e 


36 


35 42% 


D+ K^e 




37 


45 80% 



We also check consistency between isospin conjugate pairs. Isospin symmetry implies 
that total rates for K~e^Ve and K^e^Ug are approximately equal, while the 

total rate for — >■ 7r~e+fe is approximately twice that of — > n^e^Ug. After correcting 
for phase space differences, our partial rates summed over all bins agree with these 
expectations within 1.4 standard deviations. Because there are small differences in phase 
space, it is convienient to compare not rates, but form factors, as shown in Fig. [HI We obtain 
the f+{q^) at the center of bin i using 



f ( 2^ 1 Ar. 24i3 



where Ag/ is the size of bin i, \Vcd\ = 0.2256 ± 0.0010 and \Vcs\ = 0.97334 ± 0.00023 are 
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FIG. 7: Distributions of cos^vKe (upper), the cosine of the angle between the virtual W and the 
positron, positron momentum (middle) and hadron momentum (lower) in events satisfying —60 
MeV < f7 < 60 MeV. Signal and background shapes are taken from MC simulations and scaled 
using the parameters of the signal yield fits. 



from Particle Data Group fits assuming CKM unitarity 19|, and the effective in bin i 
is given by 

^ ]um'^<i! ' ^^^^ 

where f+{q^) and f+{qf) are calculated using the three parameter series parameterization 
with parameters measured in the data (see Sec. IV C\i . 

Our measured form factors in each bin are seen to be in good agreement with the 
LQCD calculations 2l|] , but with significantly smaller uncertainties, as shown in Fig. [SI 
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FIG. 8: f+{q^) comparison between isospin conjugate modes and with LQCD calculations 211 ] . 
The solid lines represent LQCD fits to the modified pole model 15|]. The inner bands show LQCD 
statistical uncertainties, and the outer bands the sum in quadrature of LQCD statistical and 
systematic uncertainties. 



The procedure for measuring partial rates is tested using the generic MC sample, from 
which events are drawn randomly to form mock data samples, each equivalent in size to the 
data sample. In each case, the measured partial rates are consistent with the input rates 
and the distributions of the deviations are consistent with Gaussian statistics. 



IV. SYSTEMATIC UNCERTAINTIES IN PARTIAL RATES 

Our determinations of the AFj are subject to a variety of systematic uncertainties. Ta- 
bles IVTIl and IVni] list each source of systematic uncertainty and its contribution to the total 
uncertainty in each of the partial rates. Because we are interested in measuring form factor 
shapes that vary with g^, it is important that we not only understand the uncertainties in 
the individual partial rates but also their correlations across q^. For each semileptonic mode 
and each significant source of systematic uncertainty, we construct an m x m (where m is 
the number of bins studied for the mode in question) covariance matrix that encapsulates 
both of these pieces of information. We now describe how each of the covariance matrices 
is estimated. 

Tag reconstruction biases enter both the numerator and denominator of our partial rate 
formulation in Eq. (jH]), and therefore largely cancel. However, there are two sources of 
systematic uncertainty related to tag yields. One source originates in the line shapes used 
to extract tag yields in data; we estimate this by using alternate line shapes and find an 
uncertainty of 0.4% for partial rates in both and D~ modes. The selection of one tag 
per mode also introduces a systematic uncertainty, primarily due to possible mismodeling 
of MC vr'' fake rates. Based on estimates of tag-fake rates in data and MC samples, we 
assign a systematic uncertainty of 0.4% to the partial rates in modes and 0.7% to those 
in D~ modes, where a greater fraction of tags contain 7r°'s. As the uncertainties associated 
with tag yields are independent of the kinematics of the semileptonic decay, they are fully 
correlated across bins. 

Systematic uncertainties associated with semileptonic track, tt", K~ , 7r°, and recon- 
struction are all studied in a similar manner: we choose fully hadronic events containing 
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TABLE VII: Summary of partial rate (AFj) uncertainties in percent for tt e'^i'e and D' 

K''e^Ve. The sign gives the direction of change relative to the change in the first bin. 





a(Ai'i) 


a(Al'2) 


\ •' 1 






o-(Al'(j) 


(t(A1'7) 


a(Ar«) 




D TT e^Ve 




















Tag line shape 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 






lag rakeH 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 






Xraekiiig efficiency 


0.48 


0.48 


0.48 


0.48 


0.48 


0.49 


0.51 






7r± ID 


0.21 


0.11 


0.05 


0.03 


0.02 


0.02 


0.04 






e± ID 


0.37 


0.38 


0.38 


0.39 


0.33 


0.18 


-0.14 






FSR 


0.18 


0.11 


0.09 


-0.02 


-0.10 


-0.20 


-0.24 






Signal shape 


0.56 


0.46 


0.58 


0.49 


0.50 


0.56 


0.49 






Backgrounds 


0.39 


0.43 


0.60 


0.61 


0.58 


0.52 


0.76 






M^C form factor 


0.06 


-0.05 


-0.05 


-0.05 


-0.07 


-0.11 


-0.04 






Q smearing 


0.84 


-0.11 


-0.26 


-0.16 


0.30 


-0.60 


-0.28 






D Lifetime 


0.37 


0.37 


0.37 


0.37 


0.37 


0.37 


0.37 






All systematic 


1.44 


1.13 


1.27 


1.22 


1.22 


1.31 


1.30 






Statistical 


6.84 


7.29 


7.90 


8.06 


8.87 


8.42 


8.63 






D" ^ K e+fe 




















Tag line shape 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


Tag fakes 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


Tracking efficiency 


0.69 


0.72 


0.75 


0.79 


0.84 


0.92 


1.04 


1.26 


1.22 


ID 


0.19 


0.13 


0.11 


0.10 


0.09 


0.10 


0.11 


0.11 


0.21 


e± ID 


0.41 


0.42 


0.43 


0.45 


0.47 


0.48 


0.44 


0.33 


0.21 


FSR 


0.12 


0.08 


0.07 


0.01 


-0.10 


-0.15 


-0.23 


-0.28 


-0.32 


Signal shape 


0.16 


0.12 


0.12 


0.14 


0.12 


0.11 


0.09 


0.14 


0.21 


Backgrounds 


0.14 


0.04 


0.12 


0.09 


0.08 


0.08 


0.04 


0.10 


0.33 


MC form factor 


0.02 


-0.02 


-0.02 


-0.01 


-0.01 


-0.01 


0.00 


0.02 


-0.08 


smearing 


0.62 


-0.11 


0.07 


-0.12 


-0.06 


-0.51 


0.08 


-0.62 


-2.05 


D Lifetime 


0.37 


0.37 


0.37 


0.37 


0.37 


0.37 


0.37 


0.37 


0.37 


All systematic 


1.26 


1.10 


1.12 


1.15 


1.20 


1.36 


1.35 


1.63 


2.55 


Statistical 


2.03 


2.19 


2.31 


2.47 


2.73 


3.14 


3.63 


4.90 


8.43 



a particle of type X, where X — 'k' ^ K~, t:^, or Kg, and reconstruct all particles in the 
event except for X. We then form missing mass squared distributions, which peak at M]^ 
for correctly reconstructed events. Wc then tally the fraction of events with the appropriate 
in which X was successfully reconstructed, after correcting for backgrounds. By doing 
this in bins of missing momentum, we compare the data and MC efficiencies as a function 
of particle momentum. 

In the case of vr and K track reconstruction and Kg finding, no evidence of bias in 
the efficiencies is found. Biases of less than 1% are observed in tt^ and K^ identification 
efficiencies. We also find 7r° reconstruction efficiencies to be approximately 6% lower in data 
than in MC simulations, the bias being roughly constant across 7r° momentum. About half 
of this discrepancy has been traced to incorrect modeling of the lateral spread of photon 
showers in the calorimeter and the energy resolution; the other half is of unknown origin. 
We reweight the MC distributions to correct for all reconstruction biases. 

Systematic uncertainties in the particle reconstruction efficiencies are often correlated 
across momentum bins. When this is the case, we construct a covariance matrix binned in 
particle momentum by noting that efficiencies binned in particle momentum (e^) are related 

2 

to efficiencies binned in semileptonic (e^^ ) via 

ei = Ae\, (13) 

where A is a matrix giving the fraction of type X particles that are part of a semileptonic 
decay in a given bin that are also in a given momentum bin; this matrix is estimated 
using signal MC simulations. The fractional covariance matrix binned in g^, , is then 
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TABLE VIII: Summary of partial rate (AFj) uncertainties in percent for — > 7r^e~''fe and 
— K^e^Ve- The sign gives the direction of change relative to the change in the first bin. 





a(Ai'i) 


a(Al'2) 








o-(Al'(j) 


(t(A1'7) 


























Tag line shape 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 






lag rakCH 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 






Tracking efficiency 


0.25 


0.25 


0.25 


0.24 


0.24 


0.24 


0.23 






tt" id 


1.06 


0.98 


1.04 


1.22 


1.83 


2.14 


1.96 






ID 


0.32 


0.32 


0.34 


0.32 


0.27 


0.13 


-0.22 






FSR 


0.14 


0.20 


0.08 


-0.05 


-0.14 


-0.22 


-0.21 






Signal shape 


1.72 


0.93 


1.91 


-1.24 


3.51 


2.43 


3.26 






Backgrounds 


0.92 


0.82 


-1.01 


0.72 


0.74 


1.38 


-6.04 






MC form factor 


0.15 


-0.03 


-0.07 


-0.06 


-0.10 


-0.15 


0.57 






Q smearing 


1.69 


0.28 


-1.74 


1.45 


-0.17 


-1.22 


-1.41 






D Lifetime 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 






All systematic 


3.01 


1.97 


3.17 


2.63 


4.18 


3.89 


-7.38 






Statistical 


9.25 


10.23 


11.24 


11.28 


13.44 


12.38 


17.98 


























Tag line shape 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


0.40 


Tag fakes 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


0.70 


Tracking efficiency 


0.76 


0.77 


0.78 


0.79 


0.81 


0.83 


0.87 


0.91 


0.96 


ID 


2.00 


1.96 


1.90 


1.83 


1.71 


1.51 


1.25 


1.35 


1.89 


e± ID 


0.42 


0.43 


0.43 


0.45 


0.48 


0.48 


0.44 


0.33 


0.20 


FSR 


0.17 


0.13 


0.08 


0.01 


-0.11 


-0.16 


-0.23 


-0.24 


-0.28 


Signal shape 


0.20 


0.22 


0.17 


0.20 


0.23 


0.26 


0.38 


0.26 


0.47 


Backgrounds 


0.13 


0.13 


0.11 


0.11 


0.14 


0.15 


0.27 


0.23 


1.46 


MC form factor 


0.03 


-0.02 


-0.02 


-0.02 


-0.02 


-0.01 


0.01 


0.02 


0.08 


smearing 


0.63 


-0.24 


-0.02 


0.29 


-1.06 


0.75 


-0.67 


-0.78 


-1.11 


D Lifetime 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


0.67 


All systematic 


2.52 


2.42 


2.36 


2.33 


2.47 


2.23 


2.08 


2.16 


3.03 


Statistical 


2.63 


2.90 


3.04 


3.23 


3.82 


3.98 


5.04 


6.88 


10.63 



given by 

= AMf A^, (14) 

where M*' is the momentum-binned fractional covariance matrix. This equation is used to 
obtain the binned systematic covariance matrices associated with track, 7r~, K~ ^ 7r°, and 
K^g reconstruction. 

The covariance matrices for all remaining systematic uncertainties - those associated with 
positron identification, FSR, background and signal shapes used to obtain signal yields, 
form factor parameterization in MC simulations, and smearing in - are estimated by 
the following procedure: for each source of systematic uncertainty, we vary the analysis in 
a manner that approximates the uncertainty on the effect in question and remeasure the 
partial rates AFj. The covariance matrix elements Afjj for this source can then be estimated 
via 

M,, = 5(Ar,)5(Ar,), (15) 

where b (AFj) denotes the difference between the partial rate in bin i measured using the 
varied analysis and the rate using the standard analysis technique. In most cases, we make 
several variations to the analysis and sum the resulting covariance matrices. Where it is 
possible to vary some parameter by positive and negative values, we average the results of 
the positive and negative variations. 

Positron identification efficiencies as a function of positron momentum are measured in 
MC simulations and in data using radiative Bhabha (ee7) and two-photon (eeee) events. 
Since the positrons in these events are relatively isolated, we embed these positrons into 
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hadronic events, and determine the decrease in efficiency. Biases of around 1.5% are ob- 
served, and the MC signal efficiency matrices and U distributions are corrected for these 
biases. To estimate the systematic uncertainty due to positron identification, we shift the 
corrections by the uncertainties on their measurement and remeasure the partial rates using 
efficiencies and U distributions obtained with the shifted corrections. 

FSR affects the partial rate measurements primarily by causing mismeasurements of 
positron momentum. FSR in the MC simulations is modeled using PHOTOS version 2.15, 
which models FSR significantly better than earher versions. To estimate the systematic 
uncertainty due to FSR, we reweight the efficiency matrices and U distributions so that 
the energy and angular distribution of photons reconstructed in the neighborhood of the 
positrons match those measured in data, and remeasure the partial rates. We have also 
studied systematic uncertainties associated with ISR, which are found to be negligible. 

The signal shapes used to model semileptonic candidates in the signal yield fits are taken 
from signal MC distributions convolved with a double Gaussian. The systematic uncertainty 
associated with this procedure is estimated by varying the widths of the Gaussians and the 
normalization of the wider one within their uncertainties and remeasuring the signal yields. 
In the D'^ — > 7r°e"*'z/e mode, there is also evidence of a possible shift between data and signal 
MC U distributions. In this mode only, we apply a systematic uncertainty equal to the 
change in rates when a shift is applied. 

The background lineshapes in the signal yield fits introduce systematic uncertainties in 
three ways. The first arises from our choice to fix the normalization of several background 
shapes, namely the small non-DD background in all modes, the p'e^Vf, and — > 

K~e^Ve backgrounds to — > n^e^v^ and the — > K^e^v^ background to Tr^e+i/g. 
The systematic uncertainties associated with these backgrounds are estimated by varying 
the normalizations within their uncertainties. In the case of the K~e^Ve and 

K^e^Ve backgrounds, where the normalizations are those that minimize the fit likelihoods 
summed over all (f' and tag modes, we vary the normalization to values that increase the 
likelihood by unity. Secondly, the choice to combine many background modes into one shape 
using fixed relative normalizations may result in incorrect background shapes. We estimate 
this systematic uncertainty by varying the normalization of several of the largest components 
of the combined shapes based on branching fraction uncertainties. Finally, incorrect MC 
fake rates may also lead to inaccurate background shapes. Our technique is most sensitive 
to positron and tt" fake rates. Using estimates of hadron-to-positron fake rates studied in 
— > i^~7r+7r+ and tt+tt" and 7r° fake rates studied in — >• i^"''7r~7r°, we estimate 

this systematic uncertainty by increasing the fake rates in MC simulations to match those 
found in data. 

The use of efficiency and smearing matrices binned in reduces the dependence of our 
results on the f+i'f) used to generate signal events in the MC simulations. However, we are 
still sensitive to non-linear effects within (f bins. To estimate the systematic uncertainty 
associated with this effect, we reweight signal MC events so that the spectra follow 
alternate form factor parameterizations. These variations sample the one standard deviation 
ellipsoid of the form factor measurements reported in this article. 

We account for possible mismodeling of the resolution in q^. Using q^ = mj^ + m\ — 
2mDEh, where Eh is the hadron energy, it follows that 5q^ = —2mDSEh. Assuming the U 
resolution is dominated by the hadron energy resolution, we estimate that the q'^ resolution 
may be 0.05 GeV^/c^ larger in data than in MC simulations in n^e^Ue and 0.02 

GeV'^/c^ larger in the other modes. Smearing the MC q^ distributions by these amounts 
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leads to the changes in the partial rates shown in Tables IVIII and IVIIII . A final systematic 
uncertainty on the partial rates arises from the and lifetimes. Using PDG 2008 
values [19i], these are 0.37% and 0.67%, respectively. 



V. FORM FACTOR, BRANCHING FRACTION, AND CKM MEASUREMENTS 

To extract form factor parameters, branching fractions, and the magnitudes of CKM 
elements \Vcd\ and \Vcs\, we fit the partial rate results using Eq. ([1]) and parameterizations 
of the form factor /_|_ (g^). Several parameterizations have been suggested. We now review 
these, before reporting the results of our fits. 



Form Factor Parameterizations 



While the exact form of /+ (g^) is not calculable in QCD, some information about the 
form factor is available. Specifically, it is expected to be an analytic fuction everywhere in 
the complex plane outside of a cut that extends along the positive axis from the mass 
of the lowest-lying cd (for D — vr) or cs (for D — > K) vector meson. This assumption leads 



to a dispersion relation (see for example Ref. |22| ) : 



u if) 



h (0) /(I - a) 



1 g 



M±(!L,u, (16) 



vr J(mr>+mp)2 t — — ie 

where rriD and mp are the masses of the semileptonic parent and daughter mesons, respec- 
tively, mD* is the mass of the D* for D ^ n ot D* for D ^ K, and a gives the relative 

contribution of this meson to /+ (0). Most of the suggested form factor parameterizations 
are motivated by a version of this dispersion relation where the integral has been replaced 
by a sum over effective poles: 



/+(0) /(!-«) 
1 - 



where the expansion parameters pk and 7^ are unknown. 

A parameterization known as the simple pole model assumes that the sum in Eq. f|T7|) is 
dominated by a single pole jl5i] : 

M<l') = (18) 
J- — — 3 — 

pole 

where the value of mpoie is predicted to be M^* . 

(s) 
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Another parameterization, known as the modified pole model 15j, adds a second term to 
the expansion given in Eq. (fTTIl . thus assuming that all higher order poles can be modeled 
by a single effective pole. To reduce the number of free parameters, this model also makes 
several simplifications, including assumptions that /?, a parameter that quantifies scaling 
violations, is near unity and 6, which describes the hard scattering of gluons, is near zero, 
leading to the prediction that 



1 + l/P-S 



m 



D 



m, 



/+(0) dq^ 



After making these simplifying assumptions, the two pole terms are reduced to 



/+(0) 



pole 



(19) 



(20) 



pole 



where mpoie is generally fixed to the D^*^-, mass and a is a free parameter. 

While the simple and modified pole parameterizations have been widely used, the presence 
of poles near the semileptonic q^ regions causes the sum in Eq. f|T7|) to have poor convergence 
properties, creating doubt as to whether truncating all but the first one or two terms leaves 
an accurate estimate of the true form factor. A third parameterization, known as the series 



expansion, attempts to address the problem [2^,123, 2J]. Exploiting the analytic properties 
of /+ (g^), a transformation of variables is made that maps the cut in the q^ plane onto a 
unit circle \z\ < 1, where 



z{q^,to) 



U~q^-Vt, 



tn 



t+-q^ + - to 



(21) 



t± = [fno ± mp)"^, and to is any real number less than t+. This transformation amounts to 
expanding the form factor about = to, with the expanded form factor given by 



1 °° r / xnfc 



(22) 



where are real coefficients, P{q^) = z{q'^,Ml),) for kaon final states, P{q^) = 1 for pion 
final states, and 0(g^,to) is any function that is analytic outside a cut in the complex 
plane that lies along the x-axis from t+ to oo. This expansion has improved convergence 
properties over Eq. f|T7j) due to the smallness of z; for example, taking the traditional 
choice oito = t+(l- {l-t-/t+y/'^y which minimizes the maximum value of z{q'^,tQ), the 
maximum value of z over the semileptonic region is 0.17 for D ^ n and 0.051 for D ^ K 



23l | . Further, taking the standard choice of 



[q^to) 



X 




-1/2 



t--q^ J 

it can be shown that the sum over all k of al is of order unity 24 . 



(23) 
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While the three parameterizations described above are the most commonly used, a fourth 
parameterization, known as ISGW2 |25|, is also occasionally used. Based on a quark model, 
this parameterization hypothesizes 



U (g^) = /. (da J (l + ^ (dax , (24) 

and predicts rigGW2 = 1-12 GeV^^c^. 



B. Fitting Technique 

Taking into account correlations across bins, our fits minimize 

m 

x'=Y: (Ar. - g{q%)Cr^\Ar, - g{q'),), (25) 

where m is the number of bins for the mode in question, Cij is the sum of the statistical 
and systematic covariance matrices for the AFj, and g{q'^)j is the predicted partial rate in 
the jth bin for the hypothesized form factor and |V^d(s)|. For each semileptonic mode, we 
perform fits using each of the four parameterizations described in Sec. IV At and provide two 
versions of the series expansion parameterization - one with only a linear term (referred to 
as the two-parameter series) and one with a linear and quadratic term (referred to as the 
three-parameter series). In all cases, we vary Vcd{s) /+(0) and one or more shape parameters: 
ri = ai/ao and r2 = 02/09 in the three-parameter series model, ri in the two-parameter 
series model, a in the modified pole model and mpoie in the simple pole model and risGW2 
in the ISGW2 model. The central values of these parameters are taken from the combined 
statistical and systematic fit. To separate the statistical and systematic uncertainties, we 
redo the fits using only statistical covariance matrices, taking the systematic uncertainty to 
be the difference between the combined and statistical-only fits in quadrature. 

If isospin is an exact symmetry, the form factors for ii^e^Ve and ir^e^Ve 

are expected to be identical, as are those for K^e^Ve and K^e^Ve- For this 

reason, we also perform combined fits to these isospin conjugate pairs. To accomplish this, 
we again minimize the given in Eq. (!25|) . now modified so that the AFj for the isospin 
conjugate pairs are combined into one vector of length 2m and Cij becomes a 2m x 2m 
covariance matrix for the combined AFj. These covariance matrices, the diagonal blocks of 
which form the covariance matrices used in the fits to individual semileptonic modes, are 
reported in the Appendix. 

The fitting technique is tested by reweighting portions of signal MC samples to alternate 
form factor parameterizations and treating these as mock data samples. We find that the 
input parameters are reproduced with no evidence of bias when the input form factor param- 
eters are the same as those used to obtain efficiency matrices. When the input parameters 
differ from the efficiency matrix paramaters, small biases are observed, and we use these to 
estimate the systematic uncertainties associated with the MC form factor parameterization, 
as described in Sec. [IVl 
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3070109-009 




2 0.5 

FIG. 9: Individual form factor fits to data (points) using 2-parameter (dashed) and 3-parameter 
(solid) series expansions. 

3070109-010 




2 0.5 

q2 (GeV2/c4) 

FIG. 10: Isospin-combined form factor fits to data (points) using 2-parameter (dashed) and 3- 
parameter (solid) series expansions. 
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TABLE IX: Results of individual form factor fits; statistical and systematic uncertainties on the 
least significant digits are shown in parentheses. For the series parameterization, we provide results 
of f+{0)Vcd, Ti = oi/oo and r2 = 02/00) as well as the expansion parameters ao, ai and 02 
themselves. The columns labeled Pij{p) give the correlation coefficients of the previous three (two) 
parameters. 



3 par. 


series 


/+(0) \Vcd(s)\ 


ri 


'■2 


Poll P02, 


P12 


X^/d.o.f. 


D° 




0.152(5)(1) 


-2.80{49){4) 


6(3){0) 


-0.44 0.68 


-0.94 


4.6/4 


D° 




0.726{8)(4) 


-2.65(34)(8) 


13{9)(1) 


-0.22 0.64 


-0.82 


3.2/6 


D+ 


■ TrOe+i/e 


0.146{7)(2) 


-1.37(88)(24) 


-4(5){1) 


-0.43 0.65 


-0.96 


0.9/4 


D+ ^ 


■ ROe+Ue 


0.707(10)(9) 


-1.66(44)(10) 


-14(11)(1) 


-0.11 0.54 


-0.82 


11.9/6 






ao 


ai 


ffl2 


POlj P02i 


P12 






n 071 (2)('[) 


-0 20(i)(()) 




-0.48 0.14 


-0.91 




DO ^ 


K-e+Ue 


n2fi4l'2V2") 


_n 07n Vo^ 




-0.22 -0.22 


-0.79 




D+ ^ 


■ TrOe+Ue 


n n74C3~n'2'i 






-0 66 37 


-0.93 




D+ ^ 


■ ROe+Ue 


0.0258(2)(3) 






-0.07 -0.26 


-0.78 




2 par. 


series 


J -T * ca{ s ) 











y2/d.o.f. 


DO ^ 




145('4~l('l') 


-1.86(18)(3) 




0.83 




8.2/5 


D" 


K-e+Ue 


7n(6)(4:) 


-2.23(19){8) 




0.68 




5.4/7 


D+ ^ 


■ nOe+u, 


0.150(5)(2) 


-1.94{25){9) 




0.80 




1.3/5 


D+ ^ 


■ ROe+Ue 


0.716{7)(9) 


-2.10{25){8) 




0.62 




13.4/7 






ao 


ffll 




P 






DO ^ 




0.071{2)(1) 


-0.13(2)(0) 




-0.89 






DO ^ 


R-e+Ue 


0.0265(2){2) 


-0.06{1)(0) 




-0.66 






D+ ^ 


■ TT^e+fe 


0.074(3)(2) 


-0.14{2)(1) 




-0.92 






D+ ^ 


• ROe+Ue 


0.0257(2)(3) 


-0.05{1)(0) 




-0.45 






Modified pole 


f+{0)\Ved(s)\ 


a 




P 




xVd.o.f. 


DO ^ 


TT e+l/e 


0.145(4)(1) 


0.21{8)(1) 




-0.82 




8.4/5 


DO ^ 


R-e+Ue 


0.716{6)(4) 


0.31{4)(2) 




-0.66 




5.9/7 


D+ ^ 


■ ttOs+u, 


0.150{5)(2) 


0.24{11){4) 




-0.77 




1.3/5 


D+ ^ 


■ ROe+v^ 


0.715{7)(9) 


0.28(6)(2) 




-0.61 




13.1/7 


Simplf 


3 pole 


f+(0)\Vea(s)\ 






P 




xVd.o.f. 


DO ^ 


TT~ e'^Ve 


0.146{3)(1) 


1.91{3)(0) 




0.71 




5.8/5 


DO ^ 


R-e+Ue 


0.720(5)(4) 


1.91{2)(1) 




0.59 




3.6/7 


D+ ^ 


■ TrOe+Ve 


0.153{4)(2) 


1.92(4)(1) 




0.65 




2.2/5 


D+ ^ 


■ ROe+Ve 


0.720{6)(9) 


1.95(3)(1) 




0.55 




14.7/7 


ISGW2 


f+{0)\Veais)\ 


r 




P 




xVd.o.f. 


DO ^ 


TT e+Ve 


0.142(4)(1) 


2.00{9)(1) 




-0.80 




12.1/5 


DO ^ 


R-e+Ue 


0.714{5)(4) 


1.60{3)(1) 




-0.64 




7.5/7 


D+ ^ 


■ TTOe+Ve 


0.148(5)(2) 


2.02(12)(5) 




-0.74 




0.9/5 


D+ 


■ ROe+Ue 


0.713{7)(9) 


1.58{4)(1) 




-0.59 




12.5/7 



C. Form Factor Results 

The optimized form factor parameters, the correlations between these parameters and the 
minimized values from fits to each semileptonic mode using each of the parameterizations 
are shown in Table lIXt the corresponding values obtained from simultaneous fits to the 
isospin conjugate pairs are shown in Table 1X1 Plots of the three parameter series expansion 
fits are shown in Figs. [9] and [TOl To facilitate display of the fit results, we have plotted 
dr/dq"^, which is estimated by dividing the AFj by the width of bin i. 

In Fig. [m we compare the form factor fits for each of the four semileptonic modes 
(ISGW2 is excluded). The partial rates (AF) obtained from each fit have been normalized 
using those from the three parameter series fits (AFseries(3))- We note that the two-parameter 
series and modified pole models give nearly indistinguishable results. Comparing the two 
and three parameter series formulations, both fits are of reasonable quality. In all modes 
but the — i> K^e~^i>e (where all parameterizations have slightly large values of due to 
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TABLE X: Results of isospin-combined form factor fits; statistical and systematic uncertainties on 
the least significant digits are shown in parentheses. For the series parameterization, we provide 
results of f-^-{0)Vcd, ri = ai/ao and r2 = 02/00) as well as the expansion parameters oq, ai and 02 
themselves. The columns labeled Pij{p) give the correlation coefficients of the previous three (two) 
parameters. 



3 par. Series 


/+(0) 




ri 


r2 


POI 


P02, 


P12 






0.150(4)(1) 


-2.35(43)(7) 


3(3){0) 


-0.43 0.67 


-0.94 


10.4/11 


D ^ K /K^e+Ue 


0.719(6){5) 


-2.25{27){7) 


3(7){1) 


-0.19 0.59 


-0.81 


19.1/15 




ao 


ai 


a2 


POl 


P02, 


P12 




D — * 'K~ I'K^e'^Ue 


0.072(2)(1) 


-0.17(3)(1) 


0.3(2)(0) 


-0.53 0.21 


-0.92 




D K~ /K°e+Ve 


0.0263{1)(2) 


-0.06(1){0) 


0.1(2)(0) 


-0.19 


-0.21 


-0.79 




2 par. Series 


/+(0) 




ri 






P 




xVd.o.f. 


D TV- /-n^e+Ve 


0.146(3)(1) 


-1.87(15)(4) 






0.81 




11.7/12 


D K-/K^>e+Ue 


0.717(4)(4) 


-2.17(15)(6) 






0.59 




19.2/16 




ao 


ai 






P 






D ■K'/TT'^e+Ue 


0.071(2)(1) 


-0.13(1){0) 






-0.89 






D K-/K°e+Ue 


0.0263{1)(2) 


-0.056(4)(2) 






-0.59 






Modified pole 


/+(0) 




a 






P 




xVd.o.f. 


D TT-/TT'>e+Ue 


0.146(3){1) 


0.21(7){2) 






-0.80 




12.2/12 


D K- /K'^e+Ve 


0.716(4)(4) 


0.30(3){1) 






-0.57 




19.4/16 


Simple pole 


/+(0) 




^polc 






P 




xVd.o.f. 


D TT~ /n^e'^Ue 


0.148(2)(1) 


1.91(2){1) 






0.68 




10.3/12 


D K-fK^e+f^ 


0.721(4)(4) 


1.93(2){1) 






0.51 




19.2/16 


ISGW2 


/+(0) 


Vcd(s) 


r 






P 




xVd.o.f. 


D 1T~ /lT^e~^Ue 


0.144(3){1) 


1.99(7){2) 






-0.78 




15.9/12 


D K-fR'^e+f^ 


0.714(4)(4) 


1.59(2)(1) 






-0.55 




20.4/16 



a statistical fluctuation between the fifth and sixth bins), the P^r degree of freedom 
using a three parameter fit is smaller than that obtained with a two parameter fit. The 
strongest evidence for a non-zero value of 02 is in n~e~^i'e, where r2 = 02/09 is slightly 

more than two standard deviations larger than zero. Thus, although there are hints of a 
preference for the three parameter fit, we do not have sufficient statistical evidence to draw 
strong conclusions on this point. 

In general, the quality of the fits is good for all parameterizations; as long as the nor- 
malization and at least one shape parameter are allowed to float, all models describe the 
data well. We take the three-parameter series fits as our nominal fits, using these to extract 
/+ (0), CKM parameters and branching fractions. We base this decision on the optimized 
convergence properties of the series expansion as well as the indications discussed above that 
the data prefer a 3-parameter fit. 

Our results rule out the predicted value of rTRGW2 and rripoie in the ISGW and simple pole 



models respectively, as have previous studies [5|, lla]. Calculating 1 + 1//5 — 5, defined in 



Eq. fll9p . using the results of the isospin-combined three parameter series fits, we find 

H- 1//3 - (5 ^ 7r"/7r°e+z/e) = 0.93 ± 0.09 ± 0.01, (26) 

and 

1 + 1//3 - 5 ^ K-/K^e+Ue) = 0.89 ± 0.04 ± 0.01. (27) 

Our values, which are in agreement with [s], [l6[ but are more precise, do not support the 
assumption by the modified pole model that 1-I-1//3 — (5~2. 

Table |XT] shows a comparison of our measurements of /-i-(0) with other experimental 
measurements and theoretical predictions. Our results are taken from the isospin-combined 
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FIG. 11: Comparison of form factor fits for each semileptonic mode. The data (squares) and fits 
to the form factor parameterizations (histograms), including the simple pole model (long dash), 
modified pole model (short dash), and two parameter series fit (dotted), are all normalized using 
the three parameter fit result (solid line at unity). 



three parameter series expansion fits and include a third uncertainty from the CKM ele- 



ments, which we take from Particle Data Group fits assuming CKM unitarity [l^. Our 
results are in good agreement with previous measurements; they are also consistent with 
LQCD predictions, although the currently available LQCD results have relatively large un- 
certainties. 

Many previous form factor measurements and predictions have used the modified pole 
model. A comparison of our measurements of the shape parameter a is shown in Table Kill 
Other experimental results are generally compatible with our results. The LQCD results 
are higher than ours by 2.1 and 2.3 standard deviations for pion and kaon final states, 
respectively. We note that variations in a between studies may reflect sensitivities to different 
regions of coupled with an imperfect parameterization rather than disagreements in the 
observed form factor distributions. Our measured form factors in each bin have been 
compared to LQCD calculations [2l|, as shown in Fig. [SI 

The BaBar experiment has reported form factor results using the series expansion for 



earlier study [5] 



They find n = -2.5 ± 0.2 ± 0.2, and ra 
16l measured ri and r2 for both K^e^u, 



0.6 ± 6.0 ± 5.0 

and 
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Our 
The 



results reported here are in agreement with previous results and are the most precise for 
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TABLE XI: Comparison of form factor normalization results to previous results. The third uncer- 
tainties show CKM uncertainties where applicable. 



/f(0) 



LQCDl [26] 
LQCD2 [21] 
Belle [27] 
BABAR [28] 
CLEO-c (281 pb- 



0.66(4)(1) 0.57(6)(2) 
0.73(3)(7) 0.64(3)(6) 

0.695(7)(22) 0.624(20)(30) 
0.727(7)(5)(7) 

[5] 0.763(7)(6)(0) 0.629(22)(7)(3) 



CLEO-c (this work) 0.739(7)(5)(0) 0.666(19)(4)(3) 



TABLE XIL Comparison of form factor shape parameter a, from fits using the modified pole 
model, with previous results. 



a 



K 



a 



LQCD [21] 
FOCUS [29] 
CLEO III [30] 
Belle [27] 
BABAR [28] 
CLEO-c (281 pb" 
CLEO-c (281 pb" 
CLEO-c (this work) 



[5] 
[16] 



0.50(4)(7) 
0.28(8)(7) 
0.36(10)(5) 
0.52(8)(6) 
0.377(23) (29) 
0.21(5)(2) 
0.21(5)(3) 
0.30(3)(1) 



0.44(4) (7) 

0.37(25)(15) 
0.10(21)(10) 

0.16(10)(5) 
0.37(8)(3) 
0.21(7)(2) 



D. Branching Fraction Results 

Branching fractions are extracted from the three parameter series expansion fit by inte- 
grating the optimized fit results over g^. We find 

7r"e+z/e) = (0.288 ± 0.008 ± 0.003)%, (28) 

^ i^-e+z/e) = (3.50 ± 0.03 ± 0.04)%, (29) 
7r°e+z/e) = (0.405 ± 0.016 ± 0.009)%, (30) 

^ /?Vz/e) = (8.83 ± 0.10 ± 0.20)%. (31) 

A comparison of these branching fractions with previous measurements is shown in Ta- 
ble IXIIII Included in the table are the averaged results of a tagged and an untagged analysis 
of the initial 281 pb^^ of CLEO-c data; differences between the results reported here and 
those of previous CLEO-c measurements are within statistical and systematic uncertainties. 
We also find that the branching fractions reported here are in excellent agreement with re- 
sults from other experiments, but are more precise. This precision arises partially from the 
ease with which D production can be determined for data collected at the ^/'(3770). 
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TABLE XIII: Comparison of branching fraction results (%) in this analysis to previous results. 



Z)0 



K-e+u, D- 



D+ ^K^ 



BES II [31] 3.82(40)(27) 

Belle \27} 0.279(27)(16) 3.45(10)(19) 

BABAR [28] 3.522(27)(45)(65) 

CLEO-c (281 pb"^) [5] 0.304(11)(5) 3.60(3)(6) 

CLEO-c (this work) 0.288(8)(3) 3.50(3)(4) 



0.378(20)(12) 
0.405(16)(9) 



8.71(38)(37) 



8.69(12)(19) 
8.83(10)(20) 



E. Extraction of \Vcd\ and \Vc< 



To extract the magnitudes of CKM matrix elements \Vcd\ and \Vcs\, we take the 



cd(s) 



/+(0) values from the isospin-combined three parameter series expansion fits and 
use the LQCD measurements [21'] /+(0) = 0.64 ± 0.03 ± 0.06 for D ^ tt transitions and 
/+(0) = 0.73 ± 0.03 ± 0.07 for D ^ K transitions. We find 



cd\ 



and 



0.234 ± 0.007 ± 0.002 ± 0.025 



0.985 ±0.009 ±0.006 ±0.103, 



(32) 



(33) 



where the third uncertainties are from /+(0). These are in agreement with those reported 
by the Particle Data Group (based on the assumption of CKM unitarity) f\^ . 



VI. CONCLUSION 

We have described measurements of the q"^ dependent partial rates of D'^ — >■ 7r~e~^i'e, 
K~e^Vei 7r°e"'"z/e, and — >■ K'^e'^Ug. We have used the partial rates to 

extract form factor parameters and branching fractions. Taking input from LQCD, we 
have reported measurements of \Vcd\ and \Vcs\- The work here uses the entire CLEO-c 
sample of ^'(3770) DD events and supersedes all previously published CLEO-c studies 
of TT'e'^i'e, — ^ K~e^Ve., ir^e^Ue-, and K^e^Ve- Our measurements 

of branching fractions and D tt~ /ir^e^Ve form factor parameters are the most precise to 
date. The results reported here are in agreement with LQCD and will be an incisive test of 
future calculations. 
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APPENDIX: CORRELATION MATRICES 



The statistical correlation matrices, described in Sec. IIIIBI and uncorrelated across the 
semileptonic modes, are shown in Tables IXIVI and [XVII The systematic correlation matrices 
are shown in Tables „ and IX VlH The diagonal blocks relating systematic uncertainties within 
a particular mode are constructed as described in Sec. IIVI To form the matrix elements for 
the off-diagonal blocks, we have assumed that the uncertainties related to tracking, tag line 
shapes, fake tags, positron identification, and FSR are fully correlated across semileptonic 
mode while all other systematic uncertainties are uncorrelated. 
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TABLE XIV: Statistical correlation matrix for — ^ tt e'^Ue and D'^ — ^ ir^e'^Ue using the standard binning, increases from left to 
right and from top to bottom. 



CO 
o 







D° TZ-e+Ue 












D+ 












1.000 -0.050 


0.001 0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 




1.000 


-0.060 0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 






1.000 -0.060 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


D° -> TT-e+I/e 




1.000 


-0.062 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 








1.000 


-0.046 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 










1.000 


-0.030 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 












1.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 














1.000 


-0.092 
1.000 


0.008 
-0.121 
1.000 


0.000 
0.011 
-0.133 


0.000 
0.000 
0.011 


0.000 
0.000 
-0.001 


-0.005 
-0.006 
-0.008 


D+ n^e+Ue 


















1.000 


-0.128 
1.000 


0.006 
-0.104 
1.000 


-0.013 
-0.016 
-0.095 
1.000 



TABLE XV: Systematic correlation matrix for — ^ tt e'^Ue and D'^ — ^ ir^e'^Ue using the standard binning, increases from left to 
right and from top to bottom. 























D+ 7r°e+i^e 










1.000 0.716 


0.595 


0.643 


0.838 


0.337 


0.409 


0.219 


0.338 


0.207 


0.241 


0.145 


0.144 


0.064 




1.000 


0.950 


0.942 


0.875 


0.831 


0.772 


0.278 


0.428 


0.264 


0.310 


0.189 


0.188 


0.083 






1.000 


0.975 


0.864 


0.906 


0.826 


0.247 


0.379 


0.235 


0.277 


0.169 


0.169 


0.074 


D° -> TT-e+I/e 






1.000 


0.912 


0.886 


0.793 


0.253 


0.384 


0.242 


0.290 


0.179 


0.181 


0.080 










1.000 


0.714 


0.709 


0.245 


0.371 


0.236 


0.286 


0.178 


0.183 


0.083 












1.000 


0.862 


0.212 


0.318 


0.205 


0.253 


0.161 


0.170 


0.083 














1.000 


0.188 


0.281 


0.181 


0.227 


0.148 


0.166 


0.092 
















1.000 


0.773 
1.000 


0.171 
0.371 
1.000 


0.614 
0.709 

-0.051 


0.675 
0.600 
0.713 


0.273 
0.625 
0.480 


-0.182 
-0.031 
0.480 


D+ n^e+Ue 




















1.000 


0.229 
1.000 


0.497 
0.536 
1.000 


0.110 
0.145 
0.273 
1.000 



TABLE XVI: Statistical correlation matrix for K e^Vg, and — > K^e^v^^ using the standard binning. increases from left to 

right and from top to bottom. 



CO 









~Ve 














D+ 












1.00 -0.05 


0.01 


0.00 0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


1.00 


-0.07 


0.01 0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 




1.00 


-0.07 0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 






1.00 -0.07 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 






1.00 


-0.06 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 








1.00 


-0.06 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 










1.00 


-0.05 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 












1.00 


-0.04 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 














1.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 
















1.00 


-0.05 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


















1.00 


-0.07 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 




















1.00 


-0.07 


0.01 


0.00 


0.00 


0.00 


0.00 






















1.00 


-0.07 


0.00 


0.00 


0.00 


0.00 
























1.00 


-0.07 


0.00 


0.00 


0.00 


























1.00 


-0.06 


0.00 


-0.00 




























1.00 


-0.06 


0.00 






























1.00 


-0.05 
































1.00 



TABLE XVII: Systematic correlation matrix for — ^ K e^v^ and — ^ K^e^v^ using the standard binning, increases from left to 
right and from top to bottom. 



CO 
CO 























D+ 




Ve 








1.00 0.80 


0.87 


0.78 0.79 


0.56 


0.81 


0.51 


0.03 


0.44 


0.46 


0.47 


0.48 


0.46 


0.51 


0.55 


0.52 


0.33 


1.00 


0.97 


0.98 0.96 


0.92 


0.92 


0.87 


0.57 


0.52 


0.54 


0.55 


0.57 


0.54 


0.60 


0.65 


0.62 


0.39 




1.00 


0.97 0.96 


0.86 


0.94 


0.82 


0.45 


0.51 


0.54 


0.55 


0.56 


0.54 


0.60 


0.65 


0.62 


0.39 






1.00 0.98 


0.94 


0.94 


0.90 


0.59 


0.51 


0.53 


0.55 


0.56 


0.54 


0.60 


0.65 


0.62 


0.39 






1.00 


0.94 


0.97 


0.90 


0.57 


0.50 


0.53 


0.54 


0.56 


0.54 


0.61 


0.66 


0.63 


0.40 








1.00 


0.89 


0.97 


0.80 


0.46 


0.48 


0.50 


0.52 


0.50 


0.56 


0.61 


0.59 


0.38 










1.00 


0.88 


0.49 


0.48 


0.51 


0.53 


0.54 


0.53 


0.60 


0.66 


0.63 


0.41 












1.00 


0.82 


0.42 


0.45 


0.46 


0.48 


0.47 


0.53 


0.59 


0.57 


0.38 














1.00 


0.26 


0.27 


0.28 


0.29 


0.29 


0.33 


0.36 


0.35 


0.24 
















1.00 


0.94 
1.00 


0.96 
0.99 
1.00 


0.98 
0.97 
0.99 
1.00 


0.75 
0.93 
0.90 
0.84 


0.96 
0.88 
0.92 
0.96 


0.77 
0.91 
0.89 
0.85 


0.75 
0.91 
0.89 
0.84 


0.66 
0.81 
0.78 
0.73 


D+ -» K°e+v^ 






















1.00 


0.70 
1.00 


0.97 
0.77 
1.00 


0.98 
0.75 
0.99 
1.00 


0.88 
0.64 
0.91 
0.91 
1.00 
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